Energy bursts in fiber bundle models of composite materials 
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As a model of composite materials, a bundle of many fibers with stochastically distributed 
breaking thresholds for the individual fibers is considered. The bundle is loaded until complete 
failure to capture the failure scenario of composite materials under external load. The fibers are 
assumed to share the load equally, and to obey Hookean elasticity right up to the breaking point. 
We determine the distribution of bursts in which an amount of energy E is released. The energy 
distribution follows asymptotically a universal power law E~^^'^, for any statistical distribution of 
fiber strengths. A similar power law dependence is found in some experimental acoustic emission 
studies of loaded composite materials. 
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I. INTRODUCTION 

During the failure process of composite materials un- 
der external load, avalanches of different magnitudes are 
produced, where an avalanche consists of simultaneous 
rupture of several elements. Such avalanches cause a sud- 
den internal stress redistribution in the material, and are 
accompanied by a rapid release of mechanical energy. A 
useful experimental technique to monitor the energy re- 
lease is to measure the acoustic emissions, the elastically 
radiated waves produced in the bursts [H, [3, 0, 11, Q ■ 

Fiber bundles with statistically distributed thresholds 
for breakdown of individual fibers are interesting models 
of failure processes in materials. They are characterized 
by simple geometry and clear-cut rules for how stress 
caused by a failed element is redistributed on the intact 
fibers. The interest of these models lies in the possibility 
of obtaining exact results, thereby providing inspiration 
and reference systems for studies of more complicated 
materials. (For reviews, see [Si,LS,i^, A The statisti- 
cal distribution of the size of avalanches in fiber bundles 
is well studied [ul, [13, [H, Q , but the distribution of the 
burst energies is not. In this article we therefore deter- 
mine the statistics of the energies released in fiber bundle 
avalanches. 

We study equal- load-sharing models, in which the load 
previously carried by a failed fiber is shared equal ly by al l 
the remaining intact fibers in the bundle (isl. 111. [l7l.[l^. 
We consider a bundle consisting of a large number N 
of elastic fibers, clamped at both ends. The fibers obey 
Hooke's law, such that the energy stored in a single fiber 
at elongation x equals ^a;^, where we for simplicity have 
set the elasticity constant equal to unity. Each fiber i is 
associated with a breakdown threshold Xi for its elonga- 
tion. When the length exceeds Xi the fiber breaks im- 
mediately, and does not contribute to the strength of 



the bundle thereafter. The individual tresholds 
assumed to be independent random variables with the 
same cumulative distribution function P{x) and a corre- 
sponding density function p{x): 



Prob(a;i < x) = P{x) 



p{v) dy. 



(1) 
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FIG. 1. The fiber bundle model. 

At an elongation x the total force on the bundle is 
X times the number of intact fibers. The average, or 
macroscopic, force is given by the expectation value of 
this, 



(F) = N X [1-P{x)]. 



(2) 



In the generic case {F) will have a single maximum Fc, a 
critical load corresponding to the maximum load the bun- 
dle can sustain before complete breakdown of the whole 
system. The maximum occurs at a critical value Xc for 
which d{F)/dx vanishes. Thus x^ satisfies 
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P{Xc) 



0. 



(3) 
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II. ENERGY STATISTICS 

Let us characterize a burst by the number n of fibers 
that fail, and by the lowest threshold value x among the 
n failed fibers. The treshold value a;niax of the strongest 
fiber in the burst can be estimated to be 



X + 



Np{x) 



(4) 



since the expected number of fibers with thresholds in an 
interval Aa; is given by the threshold distribution func- 
tion as N p{x) Ax. The last term in (0]) is of the order 
so for a very large bundle the differences in thresh- 
old values among the failed fibers in one burst are neg- 
ligible. Hence the energy released in a burst of size n 
that starts with a fiber with threshold x is given with 
sufficient accuracy as 



E 



i 71 X^ . 



(5) 



In a statistical analysis of the burst process Hemmer 
and Hansen [ll| calculated the expected number of bursts 
of size n, starting at a fiber with a threshold value in the 
interval (x, x + dx), as 



f{n, x) dx = A^- 



1 — P{x) — xp{x) 



X X{xT e^"^(^) dx, 
with the abbreviation 

X p{x) 



X{x) 



1 - P(x)' 



(6) 



(7) 



The expected number of bursts with energies less than E 
is therefore 



X! / f{n,x)dx, 



G{E) 



with a corresponding energy density 

^(2£;n)-V2/(„,^2^). 



9{E) 
Explicitly, 



dE 



g{E)=Nj2gniE), 



(8) 



(9) 



(10) 



with 



9niE) 



2E n\ 
spjs) 
1 - Pis) 



(1 - Pis) 



cxp 



sp(s)) 
l-P(s) 



(11) 



Here 



With a critical threshold value Xc, it follows from ([5]) 
that a burst energy E can only be obtained if n is suffi- 
ciently large. 



n > 2E/xl. 
Thus the sum over n starts with 

n = 1 + [2E/xl] 



(13) 



(14) 



here [a] denotes the integer part of o. 

We discuss now both the high-energy and the low- 
energy behavior of the energy density g{E). 



A. High energy asymptotics 

Bursts with high energies correspond to bursts in which 
many fibers rupture. In this range we may use Stirling's 
approximation for the factorial n!, replace H- [2E/x'j] by 
2E lx\, and replace the summation over n by an integra- 
tion. Thus 



N 



2^3/2^1/2 J „3/2 
2E/xl 



■wis) 



1 - Pis) 



exp 



(1 - Pis) - spis)) 
spis) 



1 - Pis) 



dn. 



(15) 



where s is the abbreviation lfT2|) . By changing integration 
variable from n to s we obtain 



9iE) 



N 



2^3/2^1/2 



il ^ Pis) - spis)) 



spis) ( spis) 

cxp 1 — 



1 - Pis) 
N 

2£;3/2^1/2 



ds 



I -Pis), 
(1 - Pis) - sp(s))e-^''('') 416) 



with 

his) 



1- Pis) -spis) , ^j^l--P(s) 



l-P(s) 



spis) 



\- (17) 



For large E the integral (|T6l) is dominated by the in- 
tegration range near the minimum of his). At the upper 
limit s = we have hixc) — 0, since 1—Pixc) — Xcpixc), 
Eq.|[3l). This is also a minimum of his). To see that, note 
that with y =1- sp(s)/(l - Pis)), the bracket in ifTTj) 
is of the form 



V-\nil-y)^y^+Oiy^), 



(18) 



s = y/2E/n. 



(12) 



with a minimum at y = 0. 
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In a systematic expansion about the maximum of the 
integrand in iflGj) . at s = Xc, the first factor in the integral 
(flGl) vanishes linearly, 

1 - P{s) - sp{s) = (xc ~ s)[2p(a;c) + Xcp'(xc)] 

+ 0{x,-sf, (19) 

and, as we have seen, h{s) has a quadratic minimum. 



2p{Xc) + Xcp'jXc 

xIp{xc 



(Xe - sf. (20) 



Inserting these expressions into (fTG]) and integrating, we 
obtain the following asymptotic expression. 



where 



g{E) ~ N 



C 



c 



X^p{Xc)'^ 



4^1/2 [2p{xc)+xcp'{xc)y 




(21) 



(22) 



FIG. 2. (A) The uniform threshold distribution and (B) 
the WeibuU distribution l[24|l of index 2 (dotted line) and 
index 5 (solid line) . 



In Fig. 3 we compare the theoretical formula with sim- 
ulations for the uniform distribution. 



Fix) 



X for < x < 1, 
for a; > 1, 



(23) 



which corresponds to Xc — and C — 2 ^tt and 
for the Weibull distribution with index k — 2, 

P{x) = 1 - e-""" for a: > 0, (24) 

which corresponds to Xc = 2~^/^ and C = 2^^(27re)^^/^. 
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FIG. 3. Simulation results for g{E) characterizing energy 
bursts in fiber bundles with (A) the uniform threshold 
distribution l(23|l and (B) the Weibull distribution ([24]) of 
index 2. The graphs are based on 1000 samples with A'' = 10^ 
fibers in each bundle. Open circles represent simulation data, 
and dashed lines are the theoretical results l|2H22p for the 
asymptotics. 

The corresponding asymptotics l(2T|) are also exhibited 
in Fig. 3. For both threshold distributions the agreement 
between the theoretical asymptotics and the simulation 
results is very satisfactory. The exponent —5/2 in the 
energy burst distribution is clearly universal. Note that 
the asymptotic distribution of the burst magnitudes n is 
governed by the same exponent [ll|. 
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B. Low-energy behavior 

The low-energy behavior of the burst distribution is 
by no means universal: g{E) may diverge, vanish or stay 
constant as i? ^ 0, depending on the nature of the 
threshold distribution. In Fig. 4 we exhibit simulation 
results for the low-energy part of g{E) for the uniform 
distribution and the Weibull distributions of index 2 and 
index 5. 
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FIG. 4. Simulation results for the burst distribution g{E), in 
the low-energy regime, for the uniform threshold distribution 
(circles), the Weibull distribution with k — 2 (triangles) 
and Weibull distribution with = 5 (squares). The graphs 
are based on 1000 samples with A'' = IQ^ fibers in each bundle. 

We see that g{E) approaches a finite hmit in the 
Weibull k = 2 case, approaches zero for Weibull fc = 5 
and apparently diverges in the uniform case. All this 
is easily understood, since bursts with low energy pre- 
dominently correspond to single fiber bursts (n = 1, i.e. 
E = 2:^/2) and to fibers with low threshold values. The 
number of bursts with energy less than E therefore cor- 
responds to the number of bursts with x < V^E, which 
is close to N P{V2E). This gives 



giE) c N 



V2E 



when E ^ 0. 



(25) 



For the uniform distribution g{E) should therefore di- 
verge as (2i?)~^/^ for E ^ 0. The simulation results 
in Fig. 4 are consistent with this divergence. For the 
Weibull of index 2, on the other hand, ((25l) gives g{E) 
2N when E ^ 0, a. value in agreement with simulation 
results in the figure. Note that for a Weibull distribution 
of index fc, the low-energy behavior is g{E) cx i5('=~2)/2_ 
Thus the Weibull with fc = 2 is a borderline case between 
divergence and vanishing of the low-energy density. 

The same lowest-order results can be obtained from 
the general expression ifTOj) , which also can provide more 
detailed low-energy expansions. 



III. SUMMARY 

In the present article we have studied the distribution 
of burst energies during the failure process in fiber bun- 
dles with statistically distributed thresholds for break- 
down of individual fibers. We have derived an exact ex- 
pression for the energy density distribution g{E), and 
shown that for high energies the energy density obeys a 
power law with exponent —5/2. This asymptotic behav- 
ior is universal, independent of the threshold distribution. 
A similar power law dependence is found in some experi- 
mental observations on acoustic emission studies [l|, 0| of 
loaded composite materials. 

In contrast the low-energy behavior of g{E) depends 
crucially on the distribution of the breakdown thresholds 
in the bundle. g{E) may diverge, vanish or stay constant 
for £; ^ 0. 
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